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We study the extragalactic protons with universal spectrum, which is independent of mode of propagation, 
when distance between sources is less than the propagation lengths, such as energy attenuation length or diffusion 
length (for propagation in magnetic fields). The propagation features in this spectrum, the GZK cutoff, dip and 
bump, are studied with help of modification factor, which weakly depends on the generation spectrum index 7 9 . 
We argue that from the above features the dip is the most model-independent one. For the power-law generation 
spectrum with 7 9 = 2.7 the dip is very well confirmed by the data of all existing detectors, which gives the strong 
evidence for extragalactic protons propagating through CMB. We develop the AGN model for origin of UHECR, 
which successfully explains the observed spectra up to 1 x 10 20 eV and transition from galactic to extragalactic 
cosmic rays. The calculated spectrum has the GZK cutoff, and the AGASA excess of events at E > 1 x 10 20 eV 
needs another component, e.g. from superheavy dark matter. In case of weak extragalactic magnetic fields this 
model is consistent with small-angle clustering and observed correlation with BL Lacs. 



1. Introduction 

The systematic study of Ultra High Energy 
Cosmic Rays (UHECR) started in late fifties af- 
ter construction of Volcano Ranch (USA) and 
Moscow University (USSR) arrays. During next 
50 years of research the origin of UHE particles, 
which hit the detectors, was not well understood. 
At present due to the data of the last genera- 
tion arrays, Haverah Park (UK), Yakutsk (Rus- 
sia), Akeno and AGASA (Japan), and Fly's Eye 
and HiRes (USA) we are probably very close 
to understanding the origin of UHECR, and the 
data of Auger detector [2j will undoubtedly sub- 
stantially clarify this problem. 

On the theoretical side we have an important 
clue to understanding the UHECR origin: the in- 
teraction of extragalactic protons, nuclei and pho- 
tons with CMB, which leaves the imprint on UHE 
proton spectrum, most notably in the form of the 
Greisen-Zatsepin-Kuzmin (GZK) [3] cutoff. 

We shall shortly summarize the basic experi- 
mental results and the results of the data analy- 



sis, important for understanding of UHECR ori- 
gin (for a review see 0]). 

(i) The spectra of UHECR are measured PQ with 
good accuracy at 1 - 100 EeV, and these data have 
a power to reject or confirm some models. The 
discrepancy between the AGASA and HiRes data 
at E > 100 EeV might have the statistical expla- 
nation uj. 

(ii) The mass composition at E > 1 EeV (as well 
as below) is badly known (for a review see [B]). 
The different methods give the different mass 
composition, and the same methods disagree in 
different experiments. Probably the most reli- 
able method of measuring the mass composition 
is given by elongation rate (energy dependence 
of maximum depth of shower V max ) measured by 
the fluorescent method. The data of Fly's Eye in 
1994 .1] favored iron nuclei at 1 EeV with a grad- 
ual transition to the protons at 10 EeV. The fur- 
ther development of this method by the HiRes de- 
tector, which is the extension of Fly's Eye, shows 
the transition to the proton composition already 
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at 1 EeV [7]. 

(Hi) The arrival directions of particles with en- 
ergy E > 4 x 10 19 eV show the small-angle clus- 
tering within the angular resolution of detectors. 
AGASA found 3 doublets and one triplet among 
47 detected particles jHJ (see the discussion [§]). 
In the combined data of several arrays JU| there 
were found 8 doublets and 2 triplets in 92 events. 
The stereo HiRes data do not show small- 
angle clustering for 27 events at E > 4 x 10 19 eV, 
maybe due to limited statistics. 

Small-angle clustering is most naturally ex- 
plained in case of rectilinear propagation as a 
random arrival of two (three) particles from a 
single source jT2] - This effect has been calcu- 
lated in Refs. |13ll4ll5ll6liT] . In the last four 
works the calculations have been performed by 
MC method and results agree well. According 
to ^7] the density of the sources, needed to ex- 
plain the observed number of doublets is n s = 
(1 - 3) x 10" 5 Mpc -3 . In [ini the best fit is given 
by n s ~ 1 x Mpc~ 3 and the large uncertain- 
ties (in particular due to ones in observational 
data) are emphasized. 

(iv) There have been recently found the statis- 
tically significant correlations between direction 
of particles with energies (4—8) x 10 19 eV and 
directions to AGN of the special type - BL Lacs 
[T%] (see also the criticism JH] and the reply (20j)- 
The items (Hi) and (iv) favor rectilinear prop- 
agation of primaries from the point-like extra- 
galactic sources, presumably AGN. This state- 
ment provokes many questions: 
Is rectilinear propagation the only solution to 
clustering (e.g. does it appear in propagation of 
protons in strong magnetic fields)? In case the 
correlation with AGN is true, what are the pri- 
maries (protons in weak magnetic field or neutral 
particles)? 

In this paper we shall analyse the origin of 
UHECR in two steps. In the first one we will use 
only most reliable observational data, namely the 
energy spectra and fluxes under the most conser- 
vative assumption that primaries are extragalac- 
tic protons. We shall calculate spectra for prop- 
agation in weak and reasonably strong magnetic 
fields, and demonstrate the strong evidence in fa- 



vor of protons as the primaries at energies 1 - 
80 EeV. This part of analysis is almost model in- 
dependent. In the second step we shall include 
the data (Hi) and (iv) and formulate our model. 
In the framework of our model we shall discuss 
the connection between galactic and extragalactic 
components of cosmic rays, and the problem with 
superGZK particles, i.e. ones at E > 100 EeV. 

2. Three problems of UHECR 

1. SuperGZK particles at E>lx 10 20 eV. 

"The AGASA excess", namely 11 events with en- 
ergy higher than 1 x 10 20 eV, cannot be explained 
as extragalactic protons, nuclei or photons. While 
the spectrum up to 8 x 10 19 eV is well explained 
as extragalactic protons with the GZK cutoff, the 
AGASA excess should be described as another 
component of UHECR, most probably connected 
with the new physics: superheavy dark matter, 
new signal carriers, like e.g. light stable hadron 
and strongly interacting neutrino, the Lorentz in- 
variance violation etc. 

The problem with superGZK particles is seen 
in other detectors, too. Apart from the AGASA 
events, there are five others: the golden FE event 
with E ~ 3 x 10 20 eV, one HiRes event with 
E « 1.8 x 10 20 eV and three Yakutsk events with 
E f» 1 x 10 20 eV. No sources are observed in the 
direction of these particles at the distance of or- 
der of attenuation length. The most severe prob- 
lem is for the golden FE event: with attenuation 
length Z at t = 21 Mpc and the homogeneous mag- 
netic field 1 nG on this scale, the deflection of 
particle is only 3.7°. Within this angle there are 
no remarkable sources at distance ~ 20 Mpc [JT] . 

2. Transition from galactic to extragalactic cos- 
mic rays. 

All measurements agree with the existence of the 
proton knee at energy about 2.5 x 10 15 eV and 
with increasing of the mean atomic number A of 
the primaries as energy grows up to 1 x 10 17 eV. If 
the knee is due to rigidity-dependent propagation 
(diffusion) or rigidity-dependent acceleration, the 
iron nuclei should have a knee at ~ 7 x 10 16 eV, as 
it is indeed observed by KASCADE (see Section 
6 for the details). On the other hand, the ankle 
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at E ~ 1 x 10 19 eV, discovered in late 70s in the 
Haverah Park data j^2], is traditionally consid- 
ered as transition from galactic to extragalactic 
cosmic rays. 

How the gap between 1 x 10 17 eV and 1 x 
10 19 eV is filled? 

This problem was studied in [231 - We sha11 

con- 
sider it in Section 6. 

Acceleration to E > 1 x 10 20 eV. 

Acceleration to i? mia ~ 1 X 10 21 eV is sufficient 
for present observations. It can happen that fu- 
ture observations, e.g. by EUSO , will indicate to 
considerably higher E m&x . 

The conservative shock acceleration mechanism 
can provide E max ~ 1 x I0 21 eV only for se- 
lected astrophysical objects, such as AGN [21] 
and GRBs [2S|. However, there are many other 
mechanisms of accelerations which are not devel- 
oped mathematically as good as the shock ac- 
celeration but which in principle can accelerate 
particles to much higher energies and to be very 
powerful. They include unipolar induction, oper- 
ating in the accretion discs, jets from black holes 
and fast-rotating pulsars, acceleration in strong 
e-m waves in vacuum and plasma, and different 
types of plasma acceleration mechanisms (see [2HJ 
for a review). Therefore, acceleration of particles 
to E > 1 x 10 20 eV docs not look forbidden. The 
more restrictive problem is a presence of such ac- 
celerators nearby our galaxy or inside it. How- 
ever, discovery of particles with £ > 1 x 10 20 eV 
will imply either "new" acceleration mechanisms 
or top-down scenarios. 

One can find more discussion of UHECR prob- 
lems in review |27j . 

3. Propagation theorem and the universal 
spectrum 

As numerical simulations show (see e.g. [2H1 
29 ), the propagation of UHE protons in strong 
magnetic fields changes the energy spectrum (for 
physical explanation of this effect see [SOD- The 
influence of magnetic field on spectrum depends 
on the separation of the sources d. The propaga- 
tion theorem reads: 

For uniform distribution of sources with sepa- 
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Figure 1. Convergence of the diffusive spectrum 
to the universal spectrum in the case of diffu- 
sion in the random magnetic fields with (Bq, l c ) = 
(100 nG, 1 Mpc), where Bo is the magnetic field 
on the basic turbulent scale l c . At small energies 
the Bohm diffusion is assumed. The data points 
are from AGASA. 

ration much less than characteristic lengths of 
propagation, such as attenuation length / a tt o,nd 
the diffusion length Idig, the diffuse spectrum of 
UHECR has an universal (standard) form inde- 
pendent of mode of propagation. 

For the proof see [3D] • We shall illustrate this 
theorem by convergence of the spectrum calcu- 
lated in the diffusion approximation to the uni- 
versal spectrum (see below), when the distance 
between sources, d, becomes small. In Fig. ^ the 
diffuse spectra are calculated in the diffusion ap- 
proximation for the strong magnetic field (100 nG 
on the basic scale 1 Mpc) for different distances 
between sources indicated in the figure. The uni- 
versal spectrum is shown by the solid line. All 
spectra correspond to the same emissivity C. For 
d= 50 Mpc the diffusive spectrum (black boxes) 
is quite different from the universal one. When 
d=10 Mpc (stars) both spectra are very similar, 
and at d=3 Mpc (crosses) they become indistin- 
guishable. In the case of the reasonable fields, 1 
- 10 nG, the diffusive spectra are universal for all 
reasonable separations d. 

Universal spectrum. 

One can calculate the spectrum from conservation 
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of number of particles in the comoving volume 
(protons change their energy but do not disap- 
pear). For the number of UHE protons per unit 
comoving volume, n p (E), one has: 



n p (E)dE 



dt Q g en(E g , t) dE g , 



(1) 



where t is an age of the universe, E g = E g (E,t) 
is a generation energy at age t, Q gC n(E g ,t) is the 
generation rate per unit comoving volume, which 
can be expressed through emissivity Co, the en- 
ergy release per unit time and unit of comoving 
volume at t = to, as 



Q S en(E g , t) = £ (1 + z) m Kq gen (E g ), 



(2) 



where (1 -I- z) m describes possible cosmological 
evolution of the sources. In the case of the power- 
law generation, q gcn (E g ) = E g la , with normal- 
ization constant K — "f g — 2 for 7 9 > 2. 
From Eq.QJ one obtains the diffuse flux as 
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Figure 2. Modification factor for the power-law 
generation spectra with 7 g in a range 2.0 -2.7. 
Curve r\ = 1 corresponds to adiabatic energy 
losses only, curves n ee corresponds to adiabatic 
and pair production energy losses and curves rjtot 
- to all energy losses included. 
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where 



dt/dz= H (l + z) v / n m (l + zf + n A 



(3) 



(4) 



analytic expression for dE g /dE is given in |31 |. 

The spectrum Q is referred to as universal 
spectrum. Formally it is derived from conserva- 
tion of particles and does not depend from prop- 
agation mode (see Eq. (JJJ). But in fact, the 
homogeneity of the particles, tacitly assumed in 
this derivation, implies the homogeneity of the 
sources, and thus the condition of validity of uni- 
versal spectrum is a small separation of sources. 
The homogeneous distribution of particles in case 
of homogeneous distribution of sources and inho- 
mogeneous magnetic fields follows from the Liou- 
ville theorem (see Ref. |30|). 

4. Spectrum features from proton interac- 
tion with CMB 

The extragalactic protons propagating through 
CMB produce signatures in the form of three 



spectrum features: GZK cutoff, dip and bump. 
The dip is produced due to e + e~-production and 
bump - by pile-up protons accumulated near be- 
ginning of the GZK cutoff. 

The analysis of these features is convenient to 
perform in terms of modification factor p52 33 . 

Modification factor is defined as a ratio of the 
spectrum J p (E), with all energy losses taken into 
account, to unmodified spectrum J p nTVL , where 
only adiabatic energy losses (red shift) are in- 
cluded. 



J P (E) 



7unm 
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For the power-law generation spectrum one has 



Tunm 
J p 



-in 



{ lg -2)CoE- 



7 B 



dz ^(l + z )-7 9 + l 

dz 



Modification factor is less model-dependent quan- 
tity than the spectrum. In particular, it should 
depend weakly on j g , because both numerator 
and denominator in Eq. (JjJJ include E~ lg . Fur- 
ther on we shall consider the non-evolutionary 
case m = 0. The modification factor in Fig. |2 as 
expected, depends weakly on j g , but the shape 
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Figure 3. Disappearance of bumps in diffuse 
spectra (from Ref. The sources are dis- 

tributed uniformly in the sphere of radius i?max, 
corresponding to z max . The solid and dashed 
curves are for j g = 2.7 and 7 g = 2.0, respectively. 
The curves between z max = 0.2 and z max = 2.0 
have z max — 0.3, 0.5, 1.0. 

of the GZK cutoff is strongly model-dependent: 
it is more flat in case of local overdensity of the 
sources, and more steep in case of their local 
deficit. 

The dip is a more reliable signature of interac- 
tion of protons with CMB: its shape is fixed and 
has rather complicated form to be imitated by 
other mechanism. The protons in the dip are col- 
lected from the large volume with the linear size 
about 1000 Mpc and therefore the assumption of 
uniform distribution of sources within this vol- 
ume is well justified. In contrast to this well pre- 
dicted and specifically shaped feature, the cutoff, 
if discovered, can be produced as the acceleration 
cutoff (steepening below E max ). Since the shape 
of both, GZK cutoff and acceleration cutoff, is 
model-dependent, it will be difficult to argue in 
favor of any of them. The problem of identifica- 
tion of the dip depends on the accuracy of obser- 
vational data, which should confirm the compli- 
cated shape of this feature. Do the present data 
have the needed accuracy? We shall address to 
this question in the next Section. 
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Figure 4. Predicted dip in comparison with the 
Akeno-AGASA data. 

Let us now come over to the bump. We see 
no indication of the bump in Fig. [2] at merging 
of rj ee {E) and rjtot{E) curves, where it should be 
located. The absence of the bump in the diffuse 
spectrum can be easily understood. The bump is 
clearly seen in the spectrum of a single remote 
source 32 33 . These bumps, located at different 
energies, produce a flat feature, when they are 
summed up in the diffuse spectrum. This effect 
can be illustrated by the Fig. [3| from Ref. In 
Fig. El the diffuse flux is calculated in the model 
where sources are distributed uniformly in the 
sphere of radius i? max (or z max ). When z max are 
small (between 0.01 and 0.1) the bumps are seen 
in the diffuse spectra. When radius of the sphere 
becomes larger, the bumps merge producing the 
flat feature in the spectrum. If the diffuse spec- 
trum is plotted as E 3 J p (E) this flat feature looks 
like a pseudo-bump. 

5. Dip as the signature of proton interac- 
tion with CMB 

The comparison of the calculated modifica- 
tion factor with that obtained from the Akeno- 
AGASA data, using j g = 2.7, is shown in Fig. El 
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Figure 5. Predicted dip in comparison with the 
HiRes data. 

From Fig. 01 one observes the excellent agreement 
of predicted and observed modification factors for 
the dip. By definition r/(E) < 1. In Fig. 0] one 
sees that at E < 4 x 10 17 rj b s > 1- It signals 
about appearance of another component of cos- 
mic rays, which is most probably galactic cosmic 
rays. The condition n > 1 means the dominance 
of the new (galactic) component, the transition 
occurs at higher energy. 

To calculate x 2 for the confirmation of the dip 
by Akeno-AGASA data, we choose the energy in- 
terval between 1 x 10 18 eV (which is somewhat 
arbitrary in our analysis) and 4 x 10 19 eV (the en- 
ergy of intersection of rj ee (E) and r) tot (E)). In cal- 
culations we used the Gaussian statistics for low- 
energy bins, and the Poisson statistics for the high 
energy bins of AGASA. It results in x 2 = 19.06. 
The number of Akeno-AGASA bins is 19. We 
use in calculations two free parameters: 7 S and 
the total normalization of spectrum. In effect, 
the confirmation of the dip is characterised by 
X 2 = 19.06 for d.o.f=17, or x 2 /d.o.f=1.12. 

In Fig. Othe comparison of modification factor 
with the HiRes data is shown. The agreement is 
also good. 

The good agreement of the shape of the dip 
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Figure 6. The HiRes data on the mass com- 
position. The measured X max at E > 1 x 10 18 eV 
are in a good agreement with the QGSJet-Corsika 
prediction for protons. 

T] ee (E) with observations is a strong evidence for 
extragalactic protons interacting with CMB. This 
evidence is confirmed by the HiRes data on the 
mass composition, which favor the protons (see 
Fig.\$. 

6. AGN model 

We will consider the AGN model phenomenolog- 
ically, i.e. not specifying the acceleration mecha- 
nism and assuming that space density of AGN 
satisfies the universal spectrum. The data on 
small-angle clustering and correlation with AGN 
will be involved in the analysis. We shall con- 
sider density of the sources and their luminosities, 
spectra, transition from galactic to extragalactic 
cosmic rays and the problem of superGZK parti- 
cles. 

Spectra. 

We will calculate the extragalactic proton spectra 
following rcf. 34 32], m the model with the follow- 
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Figure 7. Comparison of calculated spectra with 
the data of AGASA and HiRes. 
ing assumptions. We assume the generation spec- 
trum of a source as the standard one, oc 1/E 2 , up 
to energy E e , and more steep, oc 1/ ' E~ l9 at higher 
energies. This complex spectrum might imply 
two mechanisms of acceleration: the shock ac- 
celeration with -E max < E c and some other mech- 
anism working at higher energies. Thus, the gen- 
eration function in Eq. is determined by the 
following q gen : 



q g en(E g ) - 



E~\E g /E c 



-7s 



at Eg < E c , . 
at E g > E c , 1 > 



with normalization constant K given by 
l/K=ln{E c /E rain ) + l/{ lg -2). 



(7) 



The diffuse flux is given by Eq. (J3J. 

We consider the non-evolutionary model, m = 
0, with 7 9 = 2.7, with E c ~ 1 x 10 18 eV and 
with most conservative maximum acceleration en- 
ergy E nmx = 1 x 10 21 eV. The calculated spec- 
tra are presented in Figs. and |S] in comparison 
with the data of AGASA, Fly's Eye, HiRes and 



Figure 8. Comparison of calculated spectra with 
the data of Fly's Eye and Yakutsk. 

Yakutsk arrays. The normalization of all spec- 
tra needs somewhat different cmissivity Co- The 
normalization to the AGASA data is given by 
C = 3.5 x 10 46 erg/Mpc 3 yr. As Fig. M shows, 
the Akcno- AGASA and HiRes I - HiRes II data 
agree in the spectrum and flux, when the en- 
ergies are shifted by the factors A = 0.9 and 
A = 1.26 for the AGASA and HiRes data, respec- 
tively. Such shift is allowed by systematic errors 
in energy determination for each detector. The 
joint spectrum is htted well by the calculations 
with C = 3.1 x 10 46 erg/Mpc~ 3 yr. However, at 
E > 1 x 10 20 eV there is substantial disagree- 
ment in the data of both detectors. Statistical 
significance of this contradiction was not evalu- 
ated by the collaborations of both experiments. 
This problem will be resolved soon by the Auger 
detector. 

Transition to the galactic cosmic rays. 

In agreement with all other measurements, the 
KASCADE data (see Fig. UUJ show the gradual 
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Figure 9. Agreement between the Akeno- 
AGASA data and HiRes I - HiRcs II data 
when energies are shifted as Xagasa = 0.9 and 
^HiRes — 1-26. The combined data agree with 
calculations at E < 1 X 10 20 eV. At higher ener- 
gies the discrepancy remains. 

transition to heavy nuclei at energies above the 
proton knee E p = 2.5 x 10 6 GcV. The KASCADE 
data are in a reasonable agreement with rigidity- 
propagation models or rigidity-acceleration mod- 
els, according to which the positions of nuclei 
knees are given by Ez — ZE p , shown in Fig. ITUl 
by vertical arrows. According to this picture, 
after the iron knee, Ep c = 6.5 x 10 7 GeV, the 
CR flux should decrease steeply with energy, as 
_E~ 7 s /D(E) in the rigidity-propagation models, 
where D{E) is the diffusion coefficient in the 
galaxy. The all-particle KASCADE- Akeno spec- 
trum in Fig. EH does not show this steepening, 
and we interpret it as the compensation of the 
flux by extragalactic protons, which become the 
dominant component at E > E c ~ 1 x 10 9 GeV. 
If transition from galactic to extragalactic com- 
ponent occurs at E ~ E c ~ 1 x 10 9 GeV, it should 
reveal itself as a faint feature in all-particle spec- 
trum, because both components, galactic and ex- 
tragalactic, have similar spectra. Indeed, the dif- 
ference between spectral indices of the all-particle 
Akeno spectrum and extragalactic proton spec- 
trum of our model is equal to A7 ss 0.3. 
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Figure 10. Predicted iron nuclei spectrum (curve 
Fe model) and KASCADE data [SH]. The data 
are shown: for protons - by filled circles, for he- 
lium - by diamonds, for carbon - by stars, for 
iron - by inverted triangles, and for all-particle 
spectrum - by filled squares. The arrows labelled 
by p, He, C and Fe show the positions of corre- 
sponding knees, calculated as Ez = ZE p , with 
E p = 2.5 x 10 6 GeV. The all-particle spectrum of 
Akeno is shown by filled triangles. 



Such faint spectral feature is well known. It is 
the second knee, whose position varies from 4 x 
10 8 GeV to 8 x 10 8 GeV in different experiments 
(Fly's Eye - 4 x 10 s GeV, Akeno - 6 x 10 8 GeV, 
HiRes - 7 x 10 8 GeV and Yakutsk - 8 x 10 8 GeV). 
The second knee in the Akeno spectrum is seen 
in Fig. El 

Assuming that the total flux at E > 1 x 
10 8 GeV is given by galactic iron nuclei and ex- 
tragalactic protons, we calculated [35] the flux of 
galactic iron nuclei, subtracting the flux of ex- 
tragalactic protons, as given by our model, from 
all-particle spectrum of Akeno. The resulting 
spectra of iron nuclei are shown in Fig. El f° r 
E c = 1 x 10 9 GeV and for three values of E c in 
Fig. ^2 The latter spectra are not exactly power- 
law, but in the power-law approximation they can 
be roughly characterised by 7 rs 3.9 for spectrum 
1', 7 w 3.4 for spectrum 2', and 7 w 3.3 for spec- 
trum 3'. These spectra, especially 2' and 3' are 
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Figure 11. Calculated spectrum of extragalactic 
protons (curves 1, 2, 3) and of galactic iron spec- 
tra (curves l',2',3') compared with all-particle 
spectrum from Akeno and AGASA experiments. 
The intersections of the curves 1 — 1', 2 — 2' and 
3 — 3' give the transition from galactic (iron) to 
extragalactic (proton) components. The KAS- 
CADE data are shown by filled squares for all- 
particle fluxes and by open circles - for iron nuclei 
fluxes. 

consistent with the Hall diffusion, which predicts 
7 = 3.35 (see |3S]). 

From Fig. [*"""] one can see, that while transition 
from the galactic to extragalactic component in 
all-particle spectrum is characterised by a faint 
feature, this transition is quite sharp when the 
iron and proton spectra are resolved (the inter- 
section of curves 1' — 1, 2' — 2 and 3' — 3). 

The fraction of iron in the total flux decreases 
as energy increases from 1 x 10 17 eV to 1 x 10 18 eV, 
changing from about 80% to ~ 10% (the pre- 
dicted fraction depends on E c , see 1351). 

An interesting prediction of [!""">] is the visibil- 
ity of galactic sources at E > 1 X 10 18 eV. If 
galactic sources accelerate particles to energies 
higher than 1 x 10 18 eV, their "direct" flux can 



be seen, while the produced diffuse flux should be 
small, because of short confinement time in the 
galaxy. If generation spectrum is dominated by 
protons, the "direct" flux must be seen as the pro- 
tons, while the diffuse galactic flux is presented 
by heaviest nuclei. Due to multiple scattering 
of protons in the magnetic fields, the galactic 
point sources should be observed as the extensive 
sources with typical angular size ~ 20° at dis- 
tance r > 10 kpc. At higher energies and smaller 
distances these sources are seen as the point-like 
ones. 

SuperGZK particles from SHDM 

As Figs.[7]and[""]show, the spectrum of our model 
agrees with the HiRes, Yakutsk and Fly's Eye 
data, but does not agree with the AGASA data 
at E > 1 x 10 20 eV. The AGASA excess needs for 
its explanation another component. We will dis- 
cuss here UHECR from superheavy dark matter 
(SHDM) as a possible candidate (see also [SZj). 

SHDM as a source of UHECR without 
GZK cutoff has been suggested and studied in 
Refs. SHDM is comprised by quasi-stable 

particles with masses 10 13 — 10 14 GeV. They 
are efficiently produced at post-inflationary epoch 
(for the review see jSH])- X-particles are ac- 
cumulated in the galactic halo with overdensity 
~ 2 x 10 5 , and this effect provides the absence of 
the GZK cutoff. 

The spectra of particles produced in X-particle 
decay have been recently reliably calculated by 
three independent groups and by two different 
methods |40I41| . As the main results, these cal- 
culations give almost power-law spectrum oc E^ 1 
with 7 = 1.94 and the increased ratio of nucleons 
to photons N/j w 0.33 - 0.50 [4*1] . 

With this spectrum the SHDM model can ex- 
plain only the AGASA excess at E > 1 x 10 20 eV, 
as shown in Fig. M2L 

There are two main signatures of UHECR from 
SHDM: 

(i) The dominance of primary photons with ratio 
7 /j» w 2 - 3. 

(it) Excess of particles from the direction of the 
Galactic Center. 

The first signature has been studied recently |3""j , 
using the AGASA data. Note, that only events 
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Figure 12. The universal spectrum (dashed 
curve) and spectra from SHDM (dotted curves) 
|41| . in comparison with AGASA data. The 
SHDM spectra are shown for two normalizations. 
The sum of two components is shown by the thick 
solid curves. The x 2 values are given for the com- 
parison of these two curves with the AGASA data 
at E > 4 x 10 19 cV. 

with E > 1 x 10 20 eV are relevant. The recent 
theoretical prediction, p/-f — 0.33 — 0.50, substan- 
tially changes the conclusions of |42| . 

From 11 events at E > 1 x 10 20 eV, muons were 
detected only in G showers: 4 with two union de- 
tectors fired (cut A) and 2 more with one muon 
detector fired (cut B) (are other 5 showers muon- 
poor or muon signal in the detectors are absent 
due to fluctuations?). From Fig. 2 of Ref. 02] one 
sees only two events (cut A), where muon con- 
tent contradicts to photon-induced showers, with 
two more which are marginally consistent with 
photon-induced showers, and with two more lo- 
cated in the zone allowed for the photon-induced 
showers. With p/7 ratio given above the num- 
ber of proton-induced showers is predicted to be 
2.75 - 3.63 in a good agreement with 2 (or even 
4) showers with the hadronic muon content. As 
to arrival distribution with account of geomag- 
netic absorption and LPM suppression of show- 
ers from unabsorbed photons, we think that more 
detailed theoretical calculations and better statis- 
tics is needed for the reliable conclusions. How- 
ever, we agree with the final conclusion of the pa- 
per jl^]: "Above 10 20 eV no indication of 7-ray 



dominance is found in both /^(lOOO) and arrival 
direction distributions" . 

The second signature (ii) is reliably predicted (it 
is caused by DM distribution in the halo) and 
will be tested soon by the Auger detector. The 
absence of anisotropy centered by the Galactic 
Center excludes the SHDM model as explanation 
of the AGASA excess. 

AGN as UHECR sources. 

We will accept here the MHD simulation 03] 
of magnetic fields in various universe structures. 
These simulations favor relatively weak magnetic 
fields in the structures, namely of the order of 
0.1 nG in typical filaments and of 0.01 nG in 
the voids, the magnetic field in the Local Uni- 
verse is also weak (notice, however, the simula- 
tions 44 in which stronger magnetic fields are 
obtained) . With magnetic fields from simulations 
protons with E > 4 x 10 19 eV propagate 
quasi-rectilinearly and the small-angle clustering 
is most naturally explained with the density of 
the sources given |17I16] as 

n s ~ (1 - 3) x 10~ 5 Mpc" 3 , (8) 

with some uncertainties as indicated in ^B] . Cor- 
relation with BL Lacs [IB] indicates directly to 
AGN as the sources of UHECR. 

With n s = 3 x 10" 5 Mpc" 3 and C = 3.5 x 
10 46 erg/Mpc 3 yr, the CR-luminosity of a source 
is L p = £o/n s = 3.7 x 10 43 erg/s. Both n s and 
L p correspond well to relatively powerful AGN. 

In vicinity of Milky Way at redshift z < 0.009, 
or r < 38 Mpc, there are 12 AGN, including 
such powerful Seyferts as NGC 4051, NGC 4151, 
NGC 1068 and radiogalaxy Cen A. It corresponds 
to density n s = 5 x 10" 5 Mpc -3 , consistent 
with the density above. At redshift z < 0.0167 
(r < 70.6 Mpc), there are 19 Seyferts and radio- 
galaxies which results in n s — 1.3 x 10" 5 Mpc" 3 , 
also in agreement with the discussed density. We 
will remind that attenuation length of proton 
with E = 1 x 10 20 eV is 135 Mpc, and with 
E = 2 x 10 20 eV is 32 Mpc. The Auger detec- 
tor with large statistics at E > 1 X 10 20 eV will 
have the good chances to observe some of these 



11 



sources. 
Acceleration. 

In Eq. (jBJ we assumed the generation spectrum 
as oc l/E 2 at E < E c and cx £T 2 7 at E > E c 
with E c ~ 1 x 10 18 eV. This spectrum implies two 
mechanisms of acceleration. The shock accelera- 
tion responsible for the standard l/E 2 spectrum 
has £ max < E c . The second (hypothetical) mech- 
anism is assumed to work at E > E min ~ E c . 
As a result, the generation spectrum has some 
feature near E ~ E c , which we very approxi- 
mately describe by Eq. JJjJ. The second compo- 
nent can be due to acceleration in the jet. A 
pinch acceleration mechanism which works in jet 
plasma was suggested in Particles are ac- 

celerated with spectrum Q(E) cx E~ ls , where 
7 g = 1 + \/3 = 2.73, i.e. exactly as we assume. 
The maximum energy is connected with maxi- 
mum current at discharge E m&K — (2e/c)I max and 
can be higher than 1 x 10 21 eV. 

7. Conclusions 

We developed the most conservative scenario 
for the observed UHECR as extragalactic pro- 
tons. 

There are three signatures of UHE protons 
propagating through CMB: GZK cutoff, bump 
and dip. While bump is argued to be absent in 
the diffuse spectra and presence of the GZK cutoff 
is questioned by the AGASA data, the dip is con- 
firmed with very good accuracy by the AGASA 
and HiRes data (see Figs.^EJ). 

The predictions for the dip in terms of modi- 
fication factor are very weakly model-dependent: 
its shape varies but little with 7 g in the inter- 
val 2.0 — 2.7, the assumption of homogeneity in 
distribution of the sources is well justified, be- 
cause UHE protons are collected from the large 
distances of order l at t ~ 1000 Mpc, and its shape 
is valid for both weak and reasonably strong mag- 
netic fields. For conversion of the observed spec- 
tra into modification factor two free parameters 
are needed: j g and normalization. For 19 en- 
ergy bins of the Akeno-AGASA data and two 
free parameters, the agreement is characterised 
by x 2 = 19.06 and x 2 /d-cf.=1.12 for d.o.f.=17. 



Modification factor must satisfy rj < 1. At E < 
4 x 10 17 , as Fig. 01 shows, rj > 1, which signals 
about appearance of galactic CR component at 
energy higher, but not much than given above. 
This conclusion agrees with recent data of HiRes 
(see Fig. |SJ) which show that proton component 
becomes dominant at E > 1 x 10 18 eV. 

In our model we assume AGN as the sources, 
with generation spectrum given by Eq. JSJ. With 
emissivity Co « 3.5 x 10 46 erg/Mpc 3 yr this spec- 
trum describes well the AGASA, HiRes, Fly's Eye 
and Yakutsk data. The galactic flux of iron nu- 
clei, calculated in energy range lxlO 17 — lx 
10 18 eV from extragalactic flux of our model and 
all-particle Akeno spectrum agrees well with iron- 
nuclei flux measured by KASCADE. The spec- 
trum shape is consistent with the Hall diffusion. 

Our spectrum has the GZK cutoff and it is 
consistent with HiRes, Fly's Eye and Yakutsk 
data. For the explanation of the AGASA excess 
at E > 1 x 10 20 cV another component of UHECR 
is needed. In this paper we discuss UHECR from 
SHDM as such component. The calculated spec- 
trum agrees well with observations. We argue 
that with predicted ratio of primary photons to 
protons 7/psi 2 — 3, SHDM does not contradict 
the AGASA observations |4^]. 
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